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Quantum Monte Carlo is used to investigate the possibility of d x 2_ y 2 superconductivity in the 
two-dimensional repulsive Hubbard model. A small energy scale relevant to possible pairing requires 
a care (i.e., sufficiently small level separation between the k points (Sk, tt — 5k') and (tt — 5k" , 5k'") 
with small Sk's) to detect enhanced correlations in finite-size studies, as motivated from a previous 
study on Hubbard ladders. Our calculation indeed detects long-tailed enhancements in the d x 2_ y 2 
pairing correlation when the system is near, but not exactly at, half-filling. 



If the one-band Hubbard model, a simplest of the re- 
pulsively correlated electron systems, superconducts in 
two dimcnsions(2D), the interest is not only conceptu- 
ally generic but may be practical as well, which has in- 
deed been a challenge in the physics of high Tc supercon- 
ductivity. Some analytical calculations have suggested 
the occurrence of d^ .—2 -wave superconductivity in the 
2D Hubbard modeLErc! Numerical calculations hawe, also 
been performed extensively —Finite binding energyoEl and 
pairing interaction vertexLTO were found in those calcula- 
tions. Variational Monte Carlo calculations show thatj-a, 
superconducting order lowers the variational energyHrlLil 
Quite recently, Hotta, Takada, and the present authors 
have shown that the pairing with excluded double oc- 
cupancies has an enhanced correlationEj Nevertheless, 
there has been a reservation against the occurrence of su- 
perconductivity in the Hubbard model because the bare 
pairing correlation functions, do not show any symptom 
of long-range behavior E'Ltcj 

Now, quite a different avenue has emerged recently in 
the physics of ladders. While the. steak-coupling theory 
for the two-leg Hubbard ladderEjilJ predicts the dom- 
inant pairing correlation, the numerical calculations for 
small values of interaction gives enhancement of the pair- 
ing correlation only when the Fermi level for [/(on-site 
repulsion) = lies between a bonding and an antibond- 
ing band levels th^Lare separated with a sufficiently 
small level offset. t3~El This is because the relevant en- 
ergy scale (the spin gap in the case of ladders) is so small 
that one has to make the separation between the high- 
est occupied level (HOL) and the lowest unoccupied one 
(LUL) smaller than that to detect the pairing correlation 
in finite-size rStudies. The notion has also been success- 
fully appliedE3 to confirm the weak-coupling prediction 
that r-tha-three- leg Hubbard ladder can superconduct as 
weU.Mi3 

This view enables us to have a fresh look at the 2D 
Hubbard model. Interests here are two-fold: numerically, 
how will the QMC result behave when the care for a small 
LUL-HOL is taken. Physically, will the pair-tunneling 
mechanism remain relevant also in 2D. These are pre- 



cisely the purpose of the present study. We have indeed 
found enhancements in the bare d x 2_ y 2 pairing correla- 
tion in cases where the Fermi level lies between slightly 
separated LUL (Sk, tt — 5k') and HOL (tt — Sk", Sk), be- 
tween which the pair tunneling should occur. 

It is instructive to start with the two-leg Hubbard lad- 
der, given as, 

U= -t x ^2 ( c i,y,<r c *+i,y,<? + h - c -) 

x,y,a 

-ty ^{ c l,i,<7 C xa,a + h.c.) + U^2 n x ^n X}Vtl , 

x.a x.y 

where x(= N), y(= 1, 2), and a(—], J.) specifies the 

rungs, the chains, and the spins, respectively. According 
to the weak-coupling theory (perturbational renormaliza- 
tion + bosonization) , a gap opens in the spin excitations 
due to the relevance of interband pair tunneling between 
the Fermi points (k x , k y ) = (±k F , 0), (±fcj., tt). This con- 
comitantly makes the two-point correlation of the inter- 
chain singlet, 0^1^0^2,1 — Ci,i,lCi,2,"\ , decay slowly with 
distance. In fc-space the dominant component of this 
pair reads 

a 

where a annihilates an electron with spin a at k x = k 
in band fi(= 0, tt). 

Now, when the band structure is such that Ep inter- 
sects the bonding-band top and the antibonding-band 
bottom with k F ~ tt, k F ~ 0, intrachain nearest-neighbor 
singlet pair also has a dominant Fourier-component equal 
to eqn.(|l|) with a phase shift tt relative to the interchain 
pairing. Thus, a linear combination, cr(ci,i i<T Ci j 2,-cr — 
Ci,y,o-Ci+i,y,-a) which amounts to the d x 2_ y 2 pairing, 
should have a slow decay as well. 

Large enhancement of the inter-chain pairing correla^ 
tion has in fact been found by exact diagonalizationllS 
and by density matrix renormalization grouped when Ep 
lies close to the fc-points (0,tt) and (tt,0). Although the 
d :c 2_ !y 2-like nature of the pairing was suggested,c2l d x 2_ y 2 
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pairing correlation itself has not been calculated. So, in 
our quest for 2D, we first calculate the correlation func- 
tion with QMC. 

Here-we employ the ground-state, canonical-ensemble 
QMC,E3 where we take the free Fermi sea as the trial 
state. In most cases, we have taken the projection 
imaginary time r to be §Q/t x or larger to ensure the 
convergence. We assume periodic boundary condition 
cn+i = c\. We set t x = 1 hereafter. 

We consider the case of 56 electrons in a 30 x 2 lat- 
tice (n = 0.93) with t y = 1.975, where HOL is (0,tt) 
and LUL's are (±147r/15,0) for U — 0. We have deliber- 
ately made the LUL deviate from exactly (tt, 0) because, 
although we have stressed that the interactions across 
(tt, 0) and (0, tt) favors a d x 2_ y 2 pairing, the system be- 
comes insulating if Ef at U — (denoted by E F here- 
after) lies exactly between (tt, 0) and (0, tt), for which the 
equality k % ■ + k F = tt brings about the interband umklapp 
processes.E3Ej 

In Fig.|l|, we show the d x 2_ y 2 pairing correlation 
P(r) defined by P(r) = E^HA^r^* + Ax,y + 
Ay)0(x, y)), where 0(x, y) = J2s=±i,a v( c x, v ,<jC x+ 8, y -<j- 
Cx,v,trCx,y+S,-o), wrt h c x , 3 = c x ^,c xfi = c X)2 - P(r) for 
U = 1 has an overall enhancement over the noninteract- 
ing result, and decays slowly. 

Here we have tuned t y to make the LUL-HOL gap (de- 
noted by Ae° hereafter) as small as 0.006. This proce- 
dure is important, at least for small U , in detecting an 
enhanced pairing correlation in finite systems, ft« stressed 
in ref. ^| and in our previous publicationslij'Ej We can in 
fact see this in Fig.|], where a 2 % change into t y = 1.93 
washes out the enhancement. This may seem surpris- 
ingly sensitive, but the change in t y is accompanied by a 
more than one order of magnitude increase in Ae° —►0.1 
(without changing the U = ground state). 

Now we are in position to move on to the 2D Hubbard 
model. Here we consider the isotropic case of t y ~ t x 
and x,y — 1, •••,iV with periodic boundary condition 
in both directions, where the fc-points around (0, 7r) and 
(tt, 0) are close in energy. Our expectation from the study 
on ladders is that the pair tunneling processes between 
(5k, tt — 5k') and (tt — 5k", 5k'") may result in d x 2_ y 2 
pairing, EiJ C0S (^) — C0S (fcj/)] c kT c -kj. in 2D, but an en- 
hanced pairing correlation may be detected only when 
Ae° between those levels is small. Apart from such an ar- 
gument, the importance of the interactions around (0, tt) 
and (tt, 0) in the 2D, HubJ^d-ipodc-1 has been suggested 
by various authors. i«Jra 

In the present approach, one should be able to de- 
tect the effect of a spin (or superconducting) gap due to 
the tunneling processes in the situation where E F sit- 
uates between closely lying levels of (5k, tt — 5k') and 
(tt — 5k" , 5k'"), where Sk's are small. 
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FIG. 1. QMC result for the d x 2_ y 2 pairing correlation, 
P(r), in a 30 x 2 Hubbard ladder with 56 electrons with 
U = 1 and t y = 1.975 (□) or t y = 1.93 (O)- The error bars 
are smaller than the symbols. The dashed line represents the 
noninteracting case. 



We first take 46 electrons in 8 x 8 sites (n = 0.72) 
with t y = 0.999. For this band filling, E F lies between 
the levels of (0,±3tt/4) and (±3tt/4,0). We have taken 
t y = 0.999, because the number of electrons considered 
here would have an open shell (with a degeneracy in the 
free-electron Fermi sea) for t y = 1, which will destabilize 
QMC convergence. Taking t y — 0.999 lifts the degener- 
acy between (k\,k2) and (k2,k\) to give a tiny (< 0.01) 
but finite Ae°. 

In Fig.^| we plot the d x 2_ y 2 pairing correlation, where 
the correlation for U = 1 is clearly seen to be enhanced 
over that for {7 = especially at large distances, resulting 
in a slower decay. A similar result is shown in the inset 
for a larger system (10 x 10 with 78 electrons), where E F 
lies between (0,±47r/5) and (±47r/5,0). Also shown in 
Fig|J(a) is a result for t y = 0.95, where Ae° blows up 
to ~ 0.17, and the enhancement vanishes in accordance 
with the expectations from the ladders. 

In the above situation, LUL and HOL are taken to be 
(0, TT~5k) and (n-5k, 0) with Ae° < 0.01, while the other 
levels lie more than ~ 0.1 away from E F . One might thus 
raise a criticism that the scattering processes involving 
the states away from E F are unduly neglected. Some of 
these processes may favor the pairing, while others may 
not. We can in fact focus on the scattering processes that 
do not favor d x 2_ y 2 pairing by taking a LUL state on the 
T-M (\k x \ = \k y \) line, and situate all the others except 
a certain HOL away from E F . There, d x 2_ y 2 pairing 
should not to be favored since the d x 2_ y 2 gap function has 
nodal lines along T-M. To realize such a situation, we iru 
troduce the next nearest-neighbor hopping t' = — 0.24EZJ 
in a 10 x 10-lattice with 74 electrons where we take t y = 1 
this time.E3 E° F now lies between LUL (±2tt/5, ±2tt/5) 
along T-M line and HOL (±4tt/5, 0), with Ae° - 0.01. 
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FIG. 2. A plot similar to Fig.[j| for an 8 x 8 system with 46 
electrons, and 10 x 10 sites with 78 electrons (inset). (7 = 1 
and t y = 0.999 (□) or t y = 0.95 (CO- 



FIG. 3. A plot similar to Fig^| for an 8 x 8 system with 
60 electrons, t y = 0.999, t' = 0, and a 10 x 10 system with 
74 electrons, ty = 1, if = —0.24 (inset). In both cases (7 = 
(dashed line) or U = 1(a). 



We can see from the inset of Fig.|3J that the d x 2_ y 2 corre- 
lation is indeed no longer enhanced for U = 1. 

We have then to investigate the case where both the 
scatterings which do and do not favor d x 2_ y 2 pairing are 
taken into account on an equal footing. To realize such 
a situation, we can in fact take an 8 x 8 lattice with 60 
electrons with t y = 0.999 and t' = 0. In this case, the 
energies of the occupied levels (0,7r) and (±7r/4, ±37r/4) 
are within 0.01 to those of the unoccupied levels (tt, 0), 
(±3?r/4, ±?r/4), and (±tt/2, ±tt/2), so that the scatter- 
ings between, om., (0, 7r) and (37r/4, 7t/4), which favors 
d x 2_ y 2 pairing,E3 and those between, e.g. , (0, 7r) and 
(tt/2, 7r/2), which do not, coexist. In Fig.||(b) we can 
see that the d x 2_ y 2 correlation for U = 1 is enhanced 
over the non- interacting result, so that the coexistence 
is not detrimental. A similar fc-space configuration may 
be realized by taking 92 electrons in a 10 x 10 system 
(n = 0.92), for which an enhanced d x 2_ y 2 correlation is 
obtained as well (Fig.||). 

Our final important finding is the band-filling depen- 
dence. We have calculated S = ^2 r>3 P(r), which is a 
measure of the long-range part of the correlation, for the 
10 x 10 lattice with t y — 0.999 and t' = for various 
n = 1, 0.92, 0.78, and 0.46. Ae° is kept to be < 0.01 
throughout. The summation is restricted to r > 3 in or- 
der to eliminate short-range contributions. For n = 0.46 
(46 electrons) LUL/HOL is ((0,±3tt/5) and (±3tt/5,0)), 
which significantly deviate from (0,7r) and (ir, 0). The 
result, displayed in Fig.[|, shows that the enhancement 
in S for [7 = 1 has a maximum around a finite dop- 
ing. Although some enhancement over the [7 = result 
remains at n = 0.46, the absolute value of S becomes 
small reflecting the roundcd-off Fermi surface for small 
n. Thus the message here is that the d x 2_ y 2 pairing is 
favored near, but not exactly at, half-filling. 



To conclude, we have obtained an enhancement of 
d :c 2_ 1) 2-wave pairing correlation in the 2D Hubbard 
model. Special care has been taken on Ae° between 
the HOL (5k, n - 6k') and LUL (tt - 5k", 5k'"), as moti- 
vated from the studies on ladders. Analytical approaches 
for the 2D Hubbard model show that the relevant-en- 
ergy scale for superconductivity, As, is O(0.01i).Bla In 
this light it is natural that Ae° has to be smaller than 
~ O(0. Oli) if one wishes to detect enhanced pairing cor- 
relations in finite systems. 

Although the d x 2_ y 2 correlation is shown to be en- 
hanced only when Ae° is small, this result does not nec- 
essarily imply that a high density of states around the 
Fermi level is a prerequisite to superconductivity in the 
thermodynamic limit. Rather, we argue that the 'high 
density of states' is necessitated in finite systems to make 
the situation closer to the thermodynamic limit, where 
Ag/Ae diverges after all for any value of D(Ep)- On 
the other hand, by preferentially focusing on the scat- 
terings between, e.g., (0,7r — 5k) and (tt — 8k, 0) as in 
Fig.0, we may be mimicking the situation in many of 
the high Tq cuprates for which the photoemission stud- 
ies have revealed thajUZJ^ lies very close to the extended 
van Hove singularity!^] because the scatterings involving 
the fc-points around (0, n) and (tt, 0) would indeed be 
dominant due to the high density of states. 

Present study is restricted to relatively small values of 
U/t ~ 1. For large U/t, the strategy of paying atten- 
tion to the [7 = energy levels can fail to detect the 
enhancement of the pairing correlation, if any, because 
free-electron levels may become irrelevant. Namely, some 
energy scale Ae cS (possibly a kind of charge excitation 
energy), which vanishes in the thermodynamic limit with 
As/Ae off — > oo, should still exist in finite systems even 
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FIG. 4. The integrated pairing correlation, S, 
against band filling n for a 10 x 10 lattice with t y ■ 
U = 0(Q)orU = !(□). 
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for large U/t, but Ae cff will not be dominated by Ae°. 
Then the problem will be how to tune Ae off . 

Our results suggest that the pair-tunneling processes 
are enhanced in the 2D Hubbard model, which possibly 
leads to superconductivity. Whether such an enhance- 
ment is due to, e.g., spin fluctuations remains to be an 
open question. 

Numerical calculations were performed at the Super- 
computer Center, Institute for Solid State Physics, Uni- 
versity of Tokyo, and at the Computer Center of the Uni- 
versity of Tokyo. This work was also supported in part by 
Grant-in- Aids for Scientific Research from the Ministry 
of Education of Japan. 
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We are not arguing that t' always destroys d x 2_ y 2 pairing. 
In fact, d x 2_ y 2 pairing correlation is enhanced if E F lies 
between (5k, n — 5k') and (n — &k" , 5k'"), regardless of the 
presence/absence of t' , as will be published elsewhere. 

28 This is close to the parameter values for which the vertex 
part of the d x 2_ y 2 correlation was found to have a plateau- 
like behavior in ref. |i| Since we do not find an enhancement 
of the bare d x 2_ y 2 pairing correlation for this parameter 
set, our results should not be confused with those in ref. [], 
where the behavior of the vertex part is used as the criterion 
for superconductivity. 

29 Like in the ladders, the scatterings between exactly (it, 0) 
and (0, 7r) may not strongly favor superconductivity. We 
have numerical results suggesting such a possibility, which 
will be published elsewhere. 
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